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H. K. DrcinciQ and S. GralQ 

Bethe Center for Theoretical Physics & Physikalisches Institut der Universitdt Bonn, Nuflallee 12, 53115 Bonn, Germany 

M. K. Trenkel*! 

Max- Planck- Institut fiir Physik, Fohringer Ring 6, 80805 Miinchen, Germany 

We investigate B3 mSUGRA models, where the lightest stau, fi, is the LSP. B3 models allow for 
lepton number and R-parity violation; the LSP can thus decay. We assume one non-zero B3 coupling 
X'iji^ at Mgut, which generates further B3 couplings at Mz- We study the RGEs and give numerical 
examples. The new couplings lead to additional fi decays, providing distinct collider signatures. We 
classify the fi decays and describe their dependence on the mSUGRA parameters. We exploit our 
results for single slepton production at the LHC. As an explicit numerical example, we investigate 
single smuon production, focussing on like-sign dimuons in the final state. Also considered are final 
states with three or four muons. 
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I. INTRODUCTION 

Supersymmetry (SUSY) is one of the most 

promising extensions of the Standard Model of particle 
physics (SM) In its simplest form, we obtain the 

supersymmetric Standard Model (SSM), with a dou- 
bling of the SM particle content and one extra Higgs 
doublet. The SSM solves the hierarchy problem of the 
SM if SUSY is broken at a mass scale < 0(10 TeV). 
Therefore, SUSY should be testable at the LHC 0,11, 
which will start taking data this year. 

If they exist, supersymmetric particles are typically 
much heavier than their SM partners and at colliders 
will mostly decay rapidly. This leads to cascade decay 
chains in the detector to the lightest supersymmetric 
particle (LSP). The nature of the LSP and its pos- 
sible decay modes is thus an essential feature for all 
supersymmetric signatures. It is the purpose of this 
paper to study a novel supersymmetric phenomenol- 
ogy, namely with the lightest scalar tau (stau) fi as 
the LSP [3, |T3|. In particular we analyze in detail 
the potential ri decays in baryon-triality, B3 , models 
pU . [VX [H, [H, [il] . We then study the discovery poten- 
tial of a specific signature in this framework, namely 
resonant single slepton production at the LHC, result- 
ing in multiple muons in the final state. 



A. The B3 Framework 

The mo st g eneral renormalizable supcrpotential of the 
SSM is [II [13 
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-eatK^^Ht + -e^y^X^^kUrDfD^ . (1.3) 



Here we use the standard notation of Ref. [18| . 

The supcrpotential (|I.1|) consists of two different parts. 
Wpg, involves the lepton Y^, down-quark Yjj, and 
up-quark Yjj Yukawa matrices, which give mass to the 
leptons and quarks after elcctroweak symmetry break- 
ing. 

Wj>g , consists of lepton and baryon number violating 
operators, which tog ether can lead to rapid proton de- 
cay [H, llO, HH, Hi. The SSM thus requires an ad- 
ditional svmmetrv [ill. [T^ . [T^ to stabilize the proton. 
The most widely assumed symmetry is R-parity which 
prohibits Wy^, leading to the MSSM. But R-parity 
allows dangerous dimension-five proton decay oper- 
ators such as QQQL [2^, thus proton- hexality, Pg, 
is preferred Here, we consider a third possibil- 

ity, baryon-triality, B3. B3 is a discrete .Zs-symmetry 
which prohibits only the UDD operators in Eq. p.3|) 
but also the dangerous dimension five operators. See 
for example Refs. [H, H^, [2^ for B3 models that pro- 
vide a dark matter candidate. 

The B3 SSM has some distinguishing features com- 
pared to the MSSM [HI, H^l, which can have a strong 
impact on (hadron) collider phenomenology [H, [2^ : 

1. Lepton flavor and lepton number are violated. 

2. The renormalization group equations (RGEs) 
get additional contributions [IQ, [30|, l31|] • 

3. Neutrino masses can be generated as experimen- 
tally observed [H, [H, [sl, [H [H, [13 ■ 

4. The LSP is not stable. 
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5. Supersymmctric particles can be produced 
singly, possibly on resonance. 

Since the LSP is not stable, we are not restricted to 
the lightest neutralino xf as the LSP [3l|. The most 
general B3 SSM has more than 200 parameters and in 
principle any SUSY particle can be the LSP. Within 
the MSSM, the most widely studied constrained model 
is minimal supergravity (mSUGRA) with conserved 
Pfi and radiative electroweak symmetry breaking fsH, 
m, SH, 113, m • The 124 free parameters of the MSSM 
are reduced to only five 

A/o, A/1/2, -4o, tan/3, sgn(/^) , (1.4) 

which are fixed at the grand unification (GUT) scale, 
A/gut- We have a universal scalar mass A/q, a uni- 
versal gaugino mass M1/2, a universal trilinear scalar 
coupling Aq, the ratio of the Higgs vacuum expecta- 
tion values (vev's) tan/3, and the sign of the Higgs 
mixing parameter sgn(yLt). For a wide range of these 
parameters a Xi LSP is in fact obtained at the weak 
scale, Mz [3] . There are also wide ranges of parame- 
ter space with a fi LSP, but these are cosmologically 
excluded in the MSSM or mSUGRA 

In the B3 mSUGRA model wc consider here [lol |. 
we have six parameters at the GUT scale 

Mo, Ml/2, ^0, tan/3, sgn(/x), and A', (1.5) 

where A' stands for one non- vanishing coupling A^^j.. 
A first investigation of the parameter space has shown, 
that there are extensive regions with a neutralino, a 
stau or a sncutrino LSP [a liSl- Wc shall focus here 
on a fi LSP. fi LSP scenarios have been studied in 
the literature [§, M, M, E M, 113, Bs, 49, sl E^. 
As we now discuss, we go beyond this work in several 
aspects. 

B. New Phenomenology and Outline 

The fi LSP might decay via the dominant LiQjDk 
operator, Eq. p.3p : for example via a 4-body decay in 
the presence of a non- vanishing A'jn 

A' 

f]~ — ^ T~fi~ud. (1-6) 

An important feature of B3 mSUGRA models is that 
additional B3 couplings are generated via the RGE 
running. These new couplings can lead to 2-body de- 
cays of the fi LSP. For example, X211 will generate 
A233 which allows for the decay 

ff ^ ^I'Vr . (1.7) 

Even though A233 <C Aju, this might be the dominant 
decay mode. The decay (|L6p is suppressed by phase 
space and heavy propagators. 



We analyze in detail the conditions for a dominance 
of the 2-body decay over the 4-body decay. We pro- 
vide for the first time an extensive study of B3 fi LSP 
decays and extend and specify thus the results of [5l| , 
where a first estimate has been performed. This is 
useful when studying both pair produced and singly 
produced SUSY particles within the B3 mSUGRA 
model. Typically all heavy SUSY particle decay to 
the (fi) LSP. 

In the second half of our paper, we consider the B3 
mSUGRA model with a fi LSP and focus on resonant 
single (left-handed) charged slepton i^i and sncutrino 
i^i production at hadron colliders, which proceeds via 
a dominant LiQjDk operator: 

u,dk ^ £1,, (1.8) 

d,dk ^ M,. (1.9) 

Here, Uj (dk) is an up-type (down-type) quark of gen- 
eration j (k). 

Single slepton production allows us also to study two 
B3 couplings at a time, depending on the scenario. 
The slepton is always produced via a A' whereas the 
decay of the fi LSP in the decay chain of the slepton 
might proceed via a generated A, cf. Eq. p.7p . 

Single slepton production within a x5 LSP scenario 
leads to like-sign dileptons in the final state and has 
thus a very promi sing sig nature for experimental stud- 
ies, see Refs. [HI, l53l. Im. IHsI . [sgI ]. Here we show that 
for a fi LSP, we also obtain like-sign dilepton events 
and additionally events with three or four leptons in 
the final state. We give event rates for the LHC for 
two representative sets of B3 mSUGRA parameters. 
We also discuss the background, although a detailed 
signal over background analysis is beyond the scope 
of this paper. This is the first study of single slepton 
production in fi LSP scenarios. 

We assume in the following that only one non- 
vanishing X'^jj^ is present at AfcuT, similar to the dom- 
inant top Yukawa in the SM. Allowing for more than 
one coupling leads to stricter bounds [H, [2l|, [l^, [13, 
m, [5§|. The bounds for a single A'^-^, lie between 0(1) 
and 0(10"^) depending on the flavor indices and spar- 
ticle masses. These bounds can be up to four orders 
of magnitude stronger at Mgut if one includes the 
generation of neutrino masses [13, [3 ■ We therefore 
assume below that A'^-j, < 0(10"^) and require it to 
be consistent with the observed neutrino masses. 

Resonant slepton production at hadron colliders via 
the LiQjDk operator was first investigated in [60l.l6l|. 
using tree-level production cross sections. Three- 
Icpton final states and like-sign dilepton events were 
investigated in Ref. [H, [H, Hj, [H, ■ Ref . [bI] con- 
sidered scenarios with a gravitino LSP. Experimental 
studies by the DO collaboration at the Tcvatron were 
performed in Refs. [H, assuming a Xi LSP and a 
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THE LOW ENERGY SPECTRUM OF THE 
B3 MSUGRA MODEL WITH A n LSP 



non- vanishing X^n- The NLO QCD corrections to the 
cross section were computed in [gH, H^, [gtI [gsI . The 
SUSY-QCD corrections were included by [63|- The 
latter can modify the NLO QCD prediction by up to 
35%. In Refs. [lo, [sl, [tO, [Zl| single slepton production 
in association with a single top quark was considered. 

The outline of our paper is as follows. In Sect. |TT] 
we review the B3 mSUGRA model and approximate 
formula; for sparticlc masses. We define two B3 
mSUGRA scenarios with a fi LSP, as a refcrence for 
phcnomcnological studies. We then derive approxi- 
mate equations for the RGE generation of A from A'. 
In Sect. lIIIi we classify the different decay modes of the 
fi LSP and investigate the conditions for a dominance 
of the 2-body decay over the 4-body decay and vice 
versa. In Sect. llVi we classify all possible signatures 
for resonant single slepton production in B3 mSUGRA 
models with a fi LSP. In Sect. |V] we calculate event 
rates for like-sign dimuon events as well as for three- 
and four-muon events, at the LHC. We also discuss 
backgrounds and cuts for like-sign dimuon events. We 
conclude in Sect. I VII 



n. 



We have defined the B3 mSUGRA model in Eg^ (|L5l) 
via six input parameters at the GUT scale d, [l3| • We 
now discuss the low energy spectrum. Sparticle masses 
and couplings are obtained by running the respective 
RGEs down to the weak scale. Due to the mixing 
of different quark flavors, described by the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, the RGEs of the 
B3 couplings are not independent, but highly cou- 
pled. Therefore, a single non-zero A^^j. at the GUT 
scale generates a set of other non-zero B3 couplings 
at lower scales. Assuming a diagonal charged lepton 
Yukawa matrix Y^;, only those couplings can be gen- 
erated which violate the same lepton number as A^^j,, 
i.e. X'imn &iid Xiu. No additional source of lepton 
number violation is introduced. Phenomenologically 
particularly relevant is the generation of Xiss, which 
we discuss in detail in Sect. Ill Dl 



A. Sparticle Spectra 

The low energy SUSY particle masses depend strongly 
on the universal mSUGRAparameters p.4p and only 
weakly on A' < 0(10"^) For later use, we cite 
here approximate expressions for the relevant SUSY 
particle masses in terms of the mSUGRA parameters 
as given in [t^I , c/. also the original work in Ref . [i^l • 
The masses of the sleptons of the first and second gen- 



eration are 

m?^ = + 0.15M2/2 - sin^ OwM^ cos 2/3, 
m?^ = + 0.52M^/^ - (0.5 - sin^ 6'h')A^| cos 2/3, 
ml = + 0.52M^/2 + 0-5M| cos 2/3, 

(ILIO) 

where m^^ ^ denotes the mass of a right-/left-handed 
selectron or smuon, respectively, mp the mass of a left- 
handed electron or muon sneutrino, and 6w the elec- 
troweak mixing angle. Mz is the mass of the Z boson. 

For sfermions of the third generation, the mixing be- 
tween left- and right-handed gauge-current eigenstates 
has to be taken into account. The stau mass matrix 
squared 2H| is given by 



mi 



A 



mrBLR\ 

mrBLR ml + CRRj 



(11.11) 



with mT denoting the tau lepton mass and, expressed 
in terms of left- and right-handed third generation 
softbreaking parameters m^^ and , respectively. 



All = - (0.5 - sin^ ew)Ml cos 2/3 , 
Blr = Ar - /xtan/3 , 



(11.12) 



RR, = to|, — sin^ OwM"^ cos 2/3, 



where At is the trilinear coupling of the left- and right- 
handed stau to the Higgs. In mSUGRA, At = at 
the GUT scale. The softbreaking parameters depend 
on the mSUGRA parameters as follows [t^ . 



B3 
2 



m-f = Ml + 0.52M2/2 - -AT^, 

O 



(IL13) 



Xr = + tan^ /3) + Q.lbMl/^ + Q.iZAfj , 

where Xt parameterizes the infiuencc of the tau 
Yukawa coupling. Note, that Xt can have a strong 
impact on the stau masses due to its tan^ /3 depen- 
dence, even though Xt is suppressed by a factor 10^''. 
We will investigate this effect on the fi decay branch- 
ing ratios in the next section. 

The stau mass eigenstates fi_2 are obtained from 
the gauge eigenstates by a unitary rotation U such 
that U diagonalizes the mass matrix, UdJlfW = 
diag (m?^ , m?^ ) , yielding for the masses nif^ ^ 

1 



mi^.^ =mi + -{All + Crr) 



TyiALL-CRR.r+4mlBlj,. 



(11.14) 



The gaugino masses can be approximated in terms of 
the universal gaugino mass M1/2 [t^ . 



m^o ~ Ml = O.4IM1/2, 
m^o ~ M2 = 0.84Mi/2. 



(1L15) 



4 



Here it has been used that the hghtest neutrahno Xi 
is bino-hke in many mSUGRA models and that its 
mass can be approximated by the bino mass param- 
eter Ml at the weak scale. Accordingly, the second 
lightest neutralino X2 is mainly wino-like and its mass 
governed by the wino mass parameter M2. 

B. Reference Scenarios with a fi LSP 

For the purpose of numerical studies and as future 
reference points, wc define two specific sets of B3 
mSUGRA scenarios with a fi LSP: 

Set A: Ma = OGeV, M1/2 = 500 GeV, 

Ao = 600GeV, tan/?= 13,sgn(/i) = +1, 
a single A -^ fe ^ 0|gut, 

Set B: Ma = OGeV, A/1/2 = 700 GeV, 

Aa ^ 1150 GeV, tan /3 = 26,sgn(^) = +1, 
a single A' ,, ;, ^ 0|gut- 

(11.16) 

They a re ch osen in accordance with the following 
bounds [loH : 

• BR(B, ^ M+M") < 5.8 X 10-8 at the 95% C.L. 
obtained by the CDF collaboration [t^ . 

• 2.76 X 10-"^ < BR(6 57) < 4.34 x lO"" which 
is the central theoretical value at 2a @ using the 
experimental value of [TSj . 

• The discrepancy between experiment and the 
SM prediction of the anomalous magnetic mo- 
ment of the muon is Sa^ = a^'^P = 
(29.5 ± 8 8) X 10-^°, i.e. 3.4cr [tI, (tI Izl- 
The sets pi.l6|) are chosen such that iSa^^^^ = 
qMSSM _ gSM agrees with Ja^ within 2a. 

• Higgs mass to/jO > 112.4 GeV. This value corre- 
sponds to the LEPII bound of m/jO > 114.4 GcV 
at 95% C.L. [7^ assuming a numerical error of 
2 GcV for the mass prediction. 

• A non- vanishing coupling A-^;, at the GUT scale 
will generate a tree-level neutrino mass [l^, [12, 
HSjEIIjH^- All couplings A-^ j, in the following are 
chosen such that the tree-level neutrino mass is 
smaller than the cosmological bound on the sum 
of neutrino masses from WMAP [s^ combined 
with 2dGRFS data [H: < 0.71 eV. A 
corresponding comprehensive set of bounds for 
the mSUGRA parameter set SPSla [s^ with one 
non- vanishing coupling A^^^, is given in Ref. p^ . 
Note, that the generated trec-lcvcl neutrino mass 
depends on all mSUGRA parameters (|L5[) . The 
neutrino mass bounds on X'^^f, for Set A and Set 
B are weaker compared to those for SPSla. 





masses 


[GeV J 




masses 


[GeVJ 




bet A 


DGt 




bet A 


bet 


n 


1 '7n 

i /a 




Xi 


2U3 


zyu 


Sr 




Odd 
ZOO 


X2 


ooU 


044 


f2 


Q/in 


A p;q 

40O 


Xa 


0/1 


/ 04 


sl 


o4U 




X4 


00 / 


/DO 


Ut 




4o t 


Xi 


000 


o4y 


i>e 




4Di 


X2 


000 


/ Dl 


ii 


841 


1160 


ft" 


113 


115 


bi 


970 


1300 


H° 


643 


795 




1010 


1370 


A" 


642 


795 




1010 


1340 


H+ 


648 


799 


b2 


995 


1340 








UL 


1040 


1410 


9 


1150 


1560 



TABLE I: Sparticle masses for the B3 mSUGRA sets A and 
B as defined in Eq. HII.16p . evaluated for a renormalization 

scale Qsusy = ^ ruf^ (Qsusy) nif^ (Qsusy) using Softsusy 

2.0.10 [ill. The variation due to different X'ij^ 7^ 0|gut 
and quark mixing (see Sect. Ill C|) is below the percent level. 
The masses in the second generation coincide with those 
in the first generation. 

We use the computer programs provided by [s^, [13, [13 
to calculate BR(Bs /i+^t^), BR(6 ^ 37), and 
6a^^^ . These programs do not include the B3 cou- 
plings. But the corresponding effects are negligible for 

K,k < o{io~-') i. 

We show in Table U the supcrsymmetric mass spectra 
of the parameter sets A and B (jn.l6[) . We have ne- 
glected the mass dependence on the different non-zero 
B3 couplings which is vahd if A^^j. < 0(10-^) [§]. The 
main B3 effect on the spectrum is that we allow for a 
fi LSP. 

One naturally obtains a fi LSP spectrum for A/1/2 3> 
A/q. The large A/1/2 raises the lightest neutralino mass 
(|IL15P faster than the right-handed slepton masses 
pLlOp . It also drives the gluino and indirectly via the 
RGEs the squark masses up. We thus see in Table H] 
squark and gluino masses > 1 TeV, while the slepton 
masses are below 500 GeV. Another general feature 
of a fi LSP scenario is that the second lightest neu- 
tralino and the lightest chargino are also heavier than 
the sleptons. Therefore the only conventional super- 
symmetric decays of the left-handed sleptons are via 
the lightest neutralino. Depending on the dominant 
B3 coupling and its size, the left-handed sleptons can 
also decay into two jets. 

Nearly aU sparticles in Set B (A/1/2 = 700 GeV) 
are heavier than in Set A (A/1/2 = 500 GeV). The 
most important difference for the phenomenology at 
colliders arises from the different values of tan/3 
(tan/3 = 1 3 in Se t A, tan/3 = 26 in Set B). Accord- 
ing to Eq. (|IL13[) , the soft breaking parameters of the 
stau decrease for increasing tan /3 and thus both stau 
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mass eigenstates are reduced for large values of tan /?. 
Furthermore, the mass of the lighter stau is reduced 
due to the larger L-R- mixing, cf. Eq. pi.lip . This 
effect can be seen in Table U where the mass of the 
fi LSP is 179 GeV in Set A but only 146 GeV in 
Set B. The fi mass and tan/3 strongly influence the 
possible 2- and 4-body fi LSP branching ratios. We 
will investigate this topic in detail in Sect. IIIII 

C. Fermion Mixing 

Since the B3 RGEs are coupled, given one non-zero 
B3 coupling at Mqut, we will generate many non- 
zero couplings at the weak scale Mz ■ As we will see in 
the next section, the size of the dynamically generated 
B3 couplings depends sensitively on the composition 
of the quark Yukawa matrices. For this reason we 
prepend here a short discussion of quark mixing in B3 
models. 

Initially at Mqut, all parameters are given in the 
weak-current eigenstate basis. This includes the quark 
and lepton Yukawa coupling matrices Yu, Y^i, 
and the corresponding mass matrices niu, nid, nie. 
Since, in general, these matrices are not diagonal, we 
need to rotate the (charged) lepton and quark fields 
from the weak into the mass eigenstate basis, 

/E:F-VfL,R/ri^\ (11.17) 

with fL,R denoting the left- and right-handed fermion 
fields, respectively and VfL,R denoting the corre- 
sponding rotation matrices. The mass matrices in the 
mass eigenstate basis are then given by 

V^L niu V+j^ = diag(m„, TOc, rut), 
V^LindV+j^ = diag(md,TOs,?Ti(,), (11.18) 

VgLKle V+j^ = diag(me, 771^,771,-), 

defined at the weak scale AIz- The rotation matrices 
L R '^o* directly experimentally accessible but 
only the CKM matrix VckM) 

Vckm=V„lV+l. (11.19) 

In general, the rotation matrices for the left-handed 
fields differ from those for the right-handed fields. In 
the following, however, for simplicity and definiteness, 
we assume real and symmetric Yukawa coupling matri- 
ces, thus Vf L = Vf R. Furthermore we neglect neu- 
trino masses in this context and assume that Y^; is 
diagonal in the weak-current basis. Correspondingly, 

VeL,R = l3x3- 

To further constrain the quark Yukawa couplings, we 
restrict ourselves to the extreme cases of quark mixing 
taking place completely in the up- or the down-quark 
sector, respectively. We will refer to it as "up- type 
mixing" if 

VuL.R = VcKM, VdL^j, = l3x3, (11.20) 



at the weak scale Mz and as "down-type mixing" if 

VuL3 = l3x3, VdL.R = V+KM (11.21) 

at the weak scale. Therefore, in up- type mixing sce- 
narios, the Yukawa matrices are 

YuiMz) X v„ = ^cKM ■ diag(m„,77ic,77it) • Vckm, 

Yd{Mz) X = diag(777d, 777s,777f,), (11.22) 

and in down-type mixing scenarios, the Yukawa ma- 
trices are 

Yij{Mz) X w„ = diag(r7i„, 777^, 777*), (11.23) 
Y_d(A/z) xvd^ Vckm • diag(md, ms,mb) ■ V+p,j^, 

respectively. In the following we will consider these 
two extreme cases, (vd) is the vacuum expecta- 
tion value of the up-type (down-type) neutral CP-even 
Higgs with 

Vu = vsinf3, Vd = vcoa(3, (11.24) 

where 7; = 174 GeV is the SM vacuum expectation 
value [10^ . 

As a consequence of the non-trivial quark rotation ma- 
trices, the A'y^, couphng in Eq. jTsJ also has to be ro- 
tated from the weak basis into the quark mass basis 
for a comparison with experimental data. In case of 
up-type mixing, the LiQjDk interactions of the super- 
potential p.ip in the quark mass basis arc in terms of 
SU(2) component superfields 

A^, [N,D^; ~ E,{Y+^^),iUr]D^ . (11.25) 

In the case of down-mixing they are 

KjkWi{^GKM)]lDi'' - -Ei[/™](V(tj^j^)„fcA? . 

(11.26) 

See also Ref. [S^. However for slepton production 
cross sections, we do not take into account these 
CKM effects. If needed, the corresponding rescal- 
ing of the A' coupling can be done easily. Further- 
more the sub-dominant interactions, which include 
non-diagonal matrix elements of Vckm , do not allow 
for large production cross sections since A' enters only 
quadratically. 

D. Renormalization Group Equations 

One of the most important consequence of including 
B3 effects in SUSY models is that the LSP is no longer 
stable. This is of special interest for phcnomological 
studies if the LSP couples directly to the dominant B3 
operator. This leads to large LSP decay widths and 
to distinctive final state signatures. 

In the scenarios considered in this work, cf. Eq. (|I.5[) . 
the dominant coupling is a A^.j.; for 7 7^ 3 it does 
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not couple to the fi LSP. However, the RGEs of the 
B3 couphngs are coupled via non-diagonal entries of 
Higgs- Yukawa matrices and a A-^-j. generates dynami- 
cally other B3 couplings. Among those, we want to fo- 
cus on the Ai33 which do couple directly to the fi LSP. 



The aim of the next two sections is to study the RGEs 
of the dominant A^^j, and to quantitatively determine 
the generated Ai33. We then use these results to pre- 
dict the low energy spectrum of B3 mSUGRA scenar- 
ios given by Eq. p.5[) . We will also derive approximate 
formulae that allow for a numerical implementation of 
the running of the couplings. 

The full rcnormalization group equations for the 



B3 couplings A'yj, and Xas are 



d 



le^r'-A^fc = A^, 7i?r + A^ifc 7Qf + Kjk 1^] 



dt 

d 
'dt 



167r^-^Ai33 = Ai3/ 7!^^ + Ai;3 7^1* -t- A;33 7^* 



Li 



(11.27) 



(11.28) 



with t = In Q, Q being the rcnormalization scale. The 
anomalous dimensions 7 are listed in [loj at one-loop 
level and in [3l| at two-loop level. The RGEs simplify 
considerably under the assumption of the single B3 
coupling dominance hypothesis [13, [sO] . Products of 
two or more B3 couplings including quadratic contri- 
butions of the dominant coupling can be neglected for 
A' < 0(10^^). In this limit, the one-loop anomalous 
dimensions read 



(Yi.Y+),^. - 



3 , 8 
7:92 



7S;=2(Y+Yz, 



Li 



(11.29) 



Ei 



= 2 Y+Y 



where gi, 52, 53 are the three gauge couplings. 

From Eqs. (|II.28|) and (|II.29[) . we see that the terms 
related to 7^'^ allow for the dynamical generation of 
Ai33 by a non-zero X^^^ coupling [and vice versa for 
Eq. (11123] • AH other terms in Eq. ((TL28)) only alter 
the running of Ai33 once it is generated. The RGEs 
can be further simplified. At one-loop level, all B3 cou- 
plings but the dominant AJ-^^ and the generated Ai33 
can be neglected in the RGEs since they must be gen- 
erated first by A' and thus contribute at two-loop level 
only. 

Since we work in a diagonal charged lepton Yukawa ba- 
sis, the last term in Eq. (|II.28|) . proportional to (Y^) 




FIG. 1: Superfield diagram for the dynamical generation 
of Ai33 by at one loop order, see Eq. (|II.32p . 



does not contribute to the running of Ai33. It is only 
non-zero if z = 3, but owing to the zj-antisymmetry of 
Xijk no coupling is generated in this case (A333 = 0). 

Next, a general orderin g of the parameters in the 
anomalous dimensions is [l03l | 



53 



■2> 



> <?2 > 5? 



(11.30) 



and all other entries of the Y matr ices are smaller 
by at least one order of magnitude |l04l |. The con- 
tributions to the RGEs are thus largest for diagonal 
anomalous dimensions. 

As a result, the RGEs for a non-zero X'^j^. at the GUT 
scale and a generated Ai33 reduce to 



= A' 



ijk 



7 2 o 2 16 2 
Y^5i - 3.92 - y33 



(Yl.)33(243 



3j3 



Yt/)33^.3 + (Y^)33<5. 



16- ^A.33 



= A. 



i33 



3A' 



ijk 



33"J3 

9 

5 



i3 



33 



(n.3i) 

(11.32) 



A similar analytical approximation for the generation 
of A is derived in [5l[. But the effect of the gauge 
couphngs is neglected there. See also Rcf. jsij . 

The last term in Eq. (|II.32|) induces the dynamical 
generation of Ai33. Diagrammatically, this process 
can be understood as shown in Fig. [TJ We see that 
at one-loop the lepton-doublet superfield mixes with 
the Higgs doublet superfield Hd via the B3 coupling 
X'^jf. and the standard down quark Yukawa coupling. 

Hd then couples via the tan Yukawa coupling (Yb)^^ 
purely leptonically. The resulting effective interaction 
is of the Ai33-typc. 

It is important to notice that the generation is related 

^i33 



to Yr 



Ijk- 



Whether a given X[^f. can generate A 



or not depends on whether {Yu) ^ 0. For j ^ k it 
thus depends crucially on the origin of the CKM mix- 
ing: is it dominantly down-type or up-type mixing. In 
case of down- type mixing, all entries of the Y^ matrix 
are non-zero and all A^^ j. can therefore generate a Ai33. 
In contrast, if the quark mixing takes place in the up- 
sector, only the diagonal entries of Y^ arc non-zero 
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and j = A; is required. The flavor and size of the gen- 
erated coupling depends on tan/3 and on the precise 
j, k configuration. A strong ordering is expected that 
goes along with the ordering of the entries of the Yjj 
matrix. 

In order to study the running of the B3 couplings, 
the RGEs for the Yukawa matrix elements (Yc)^.^, 

(Yjj)^^, and (Ye)^^ and the gauge couplings are also 
needed. The full RGEs for the Yukawa couplings are 
given in [l^, HOl ■ Applying the single coupling domi- 
nance hypothesis, neglecting quadratic terms in A^^j,, 
and considering only the dominant terms Eq. (jll.30p . 
they read 



+ 6{Yu)l + {Yn)l, 



16 2 

^53 



16-'|(Y^.)33=(Y^)33 



167r' 



9 2 



(Yz,)^.J-^5?-3g| 



(11.33) 



(11.34) 

16 2 
y53 



- 243) 



^{Yu)Is,, + {Ye)1]. (11.35) 



The one-loop order RGEs for the three gauge couplings 
within the MSSM arc given by ^ 




non-zero coupling at the GUT scale: 

?i'(233) = 0.01 

l'(232,231) = 0.01 

?i'(21 3,223) = 0.01 
r(21 1 ,212,221 ,222) = 0.01 



10' 



,10- 



a> 



c 



10" 



10 ' 



10' 



10' 



down-type mixing 




non-zero coupling at the GUT scale: 

E —X'{233) X'(222) 

--A,'(232) X,'(223) X'(21 2,221) 

-- X,'(231) X'(213) X'(211) 



10-^ 



10"^ 



10° 10" 
scale Q [GeV] 



10^ 



10' 



FIG. 2: Running of B3 couplings assuming a single non- 
zero A' = 0.01 coupling at the GUT scale (upper panel) 
leading to a non-zero A233 coupling (lower panel) at lower 
scales within the B3 mSUGRA scenario Set A for down- 
type mixing. 



-hgf 



(11.36) 



with bi = {33/5, 1, -3} for i = 1, 2, 3. Thus in total 
a set of nine coupled differe ntia l equations, Eqs. pi.31|) 
- pL36)) . has to be solved [lOst . 



E. Numerical Results 

For the numerical implementation of the RGEs we 
start from the framework provided by Softsusy 
2.0.10 [S^. First, Softsusy evaluates all necessary 
parameters at the SUSY scale 

Qsnsy = \J (Qsusy) (Qsusy) ■ (11.37) 

In a second step, we apply the (R-parity conserving) 
RGEs (jII.33P - (jII.36[) to run the Yukawa couplings and 
gauge couplings up to the GUT scale. Here we add 
the B3 couphngs A-^^. ^ 0|gut and Ai33 = 0|gut and 
evolve these couplings down to the scale Q using the 
above given B3 RGEs ((TL3T1) and (|IL32|) . We have 
implemented the RGEs using a standard Rungc Kutta 
formalism [Olj . 

In Figs. [21 [3l we show the running of different A2 ^, cou- 



plings, starting with A^ 



ijk 



0.01 



GUT: 



2jk 

for the case of 



down- and up-mixing respectively. In the correspond- 
ing lower panel, we show the scale dependence of the 
generated A323 = — A233 coupling. Here, we use the 
mSUGRA parameters of Set A (tan/3 = 13). 

We see that the dominant A-^j. coupling grows by 
about a factor of 3, running from the GUT scale to 
the weak scale. This effect is mainly due to the gauge 
couplings, see Ref. [Slj, where the Yukawa couplings 
were omitted. Including the Yukawa couplings reduces 
this effect, maximally for j = fc = 3. The gener- 
ated A233 coupling is at least two orders of magnitude 
smaller than the original A' coupling. Furthermore 
it depends sensitively on the flavor structure [ijk) of 
the original A' coupling. This reflects the dependence 
on the Yukawa matrix (Y/))^.^,. In case of down-type 
mixing, the ordering of the corresponding entries is 



Yl.)33 > (Yi 



'23,32 



'22 



> 



> (Yd) 12,21 > (Yi?) 13,31 > (Yd) IP 

(n.38) 

reflecting precisely the ordering of the generated cou- 
plings in Fig. [21 Small differences between the cou- 
plings generated by A^23 iKiz) or K-i2 (^i3i) are re- 
lated to the different running of the respective A' and 
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_g _ non-zero coupling at the GUT scale: 

E V(233) 

X\222) 

V(211) 

9 U lA lA lA lA ij ij lJ lA lA 



10 



\ 



10^ 



10 



10° 10" 
scale Q [GeV] 



10^* 10^'' 



F. Comparison with the Program Softsusy 

In this section, we compare our results for A'^-^, and the 
generated couphng A^sa at the SUSY scale, Eg. (|II.37p , 
with an unpublished version of Softsusy j92 |. This 
version of Softsusy contains the complete one loop 
RGBs for A^^fc (|Ili7| and A^gs pL28)l . without our 
approximations. 

We show in Table |TT] our results and the results of 
Softsusy for the case of down-type mixing and up- 
type mixing assuming different couplings A^^^. =0.01 
at the GUT scale. For the other parameters, we con- 
sider the Set A of Eq. ((TLT6)) . 

At the SUSY scale, the differences between our re- 
sults and Softsusy for the case of down-type mix- 
ing, are less than 2% for all A'^-j, couplings and less 
than 4% for the A^aa, respectively. In case of up-type 
mixing, we find the same for the couplings A^^j. with 
j = k. However for j ^ k and up-type mixing, we ob- 
serve a discrepancy between our results and Softsusy 
for the coupling A233 generated by A223 7^ 0|gut and 
-^231 7^ 0|gut, respectively. This behavior can easily 
be understood. 



FIG. 3: Same as Fig. [J] but for quark mixing in the up- 
sector. 

{Y d) ^1. coupling, depending in turn on whether j or 
k equals 3. 

In the case of up-type mixing. Fig. [31 not all A' cou- 
plings can generate a A. Since the down Yukawa cou- 
pling is diagonal, j = k is required. Other couplings 
can generate Ai33 at higher loop levels only and are 
not included in our approximations. 

Our results can easily be translated to other scenar- 
ios: The running of the dominant coupling A' is mainly 
driven by gauge interactions, Eq. pi.3ip . and thus de- 
pends only weakly on the specific SUSY parameters. 
The dependence of the generated coupling A on SUSY 
parameters is more involved but we expect tan/3 to 
have the largest impact. In general, the generated 
A coupling scales with tan^ /3, 

A,33 (X tan^ /3 , (11.39) 

if tan^ (3^1. This is because the down-quark Yukawa 
couplings {Y [and the tan Yukawa coupling 

(Ye)^^ are proportional to l/cos/3 = ^1 + tan^ /3, 
which directly follows from Eqs. pi.22p - pi.24|) . There- 
fore the magnitude of the generated A coupling for 
other scenarios can be estimated by rescaling A of 
Fig. [5] and Fig. [3] according to Eq. firSQ)) . 



The off-diagonal Yukawa matrix elements (YD)jk are 
equal to zero at the weak scale for up-type mixing. 
Running from the weak scale to the GUT scale gener- 
ates Yukawa couplings (Yr))jk, j 7^ fc, at the one loop 
level [l3,[33|. The generation of A233 via Eq. (|lL32|l oc- 
curs therefore formally at two-loop level and has been 
neglected in our approximation. In Softsusy this two- 
loop effect is taken into account and small couplings 
are generated also for j ^ k and up- type mixing. 
Compared to the case of down-type mixing, see Ta- 
ble [III the A233 couplings are suppressed by five (with 
•^231 = 0.01|gut) and three (with A223 = 0.01|gut) or- 
ders of magnitude. Note that the generation of {Y D)jk 
is not the only two loop effect that enters the full RGEs 

[13, Him. 

Therefore, our approximation for the generation of 
Ai33 by a non-zero A-^^, at the GUT scale pi.32p breaks 
down in the case of up- type mixing and j ^ k. But 
concerning fi LSP decays, the corresponding 2-body 
decay branching ratio for Ai33 is negligible compared 
to the 4-body decay branching ratio via A-^^. and our 
approximations are applicable for such phenomenolog- 
ical studies. For example, the 2-body decay branch- 
ing ratio for up-type mixing and A231 = 0.01|gut or 
A223 — 0.01 1 GUT is less than 10~^ in Set A. 

We conclude that our approximations are valid for the 
signal and decay rates that we study in this work. We 
also note that we have provided an independent check 
of the yet-to-be published version of Softsusy (o^ . 
Using a different set of mSUGRA parameters leads to 
a similar level of agreement. 
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TABLE II: Comparison between our results, Eq. (|II.3ip and Eq. (|II.32|l . and the results of an unpublished version of 
Softsusy [s^l for \[jf, and the generated couphng Ai33 at the SUSY scale, Eq. (|II.37p . We choose different couplings 
Kjk = 0.01 at the GUT scale as given in the first column of the table. The running of X[j^. is the same for down- and 
up-type quark mixing. The generation of Ai33 depends on the quark mixing assumptions and the values at the SUSY 
scale are given separately. The remaining mSUGRA parameters are these of Set A (|11.16|l . 



III. fi LSP DECAYS IN B3 MSUGRA 



A. General LSP Decay Modes 

As WG showed in Sect. [Ill a non- vanishing coupling A^^^ 
at the GUT scale generates an additional coupling Ai33 
at the weak scale which is roughly at least two orders of 
magnitude smaller than A'y J,, cf. Figs. [2] and [3l In this 
section, we compare the possible decay modes of the 
LSP via these two couplings for different B3 scenarios. 

First, let us discuss Xi LSP scenarios. The leading 
order decay modes of the Xi LSP via the dominant A 
and the generated Aiaa couplings are all three body 
decays, 



_0 ^'jk 

Xi — ' 



if Uj dk 




'"'^'^ , (III.40) 



and 



~0 -'^i33 

Xi — ' 



~0 -^is; 

Xi — 



(III.41) 



The corresponding partial widths depend quadrati- 
cally on A'^^^ and Aiss, respectively [sl, [ol, [H, [96l |. 
Therefore, the Xi decay via Ai33 is heavily suppressed 
and a Xi LSP decays predominantly via A'-^^ into SM 
particles. 

The situation changes if one considers B3 mSUGRA 
scenarios with a fi LSP, where the fi couples not di- 



rectly to the LiQjDk operator, i.e. 



1,2. In this 



case, the ri must first couple to a virtual gaugino. 
The gaugino then couples to a virtual sfermion which 
then decays via A'^^., resulting in a 4-body decay of 
the fi LSP. The possible decay modes via a virtual 
neutralino are 



T~ if Uj dk 
T~ i~ Uj dk 

dj dk 
Vi dj dk 



(III.42) 



4-body decays via a virtual chargino are also possible 
but they are suppressed due to the higher chargino 
mass in comparison to the lightest neutralino mass, 
'^(xf) > "^(Xi)- Furthermore, the (mainly right- 
handed) fi LSP couples stronger to the (bino-like) 
lightest neutralino than to the (wino-like) lightest 
chargino. 

On the other hand, the ti can directly decay via Ai33 
into only two SM particles 



T~ Vi 
17 V. 



(III.43) 



We show in Fig. [?] (Fig. O, example diagrams for the 
4-body (2-body) decay of a f\ LSP via Aj^^, (A233). 
Although the 2-body decay suffers from the small cou- 
pling, the 4-body decay is phase space suppressed as 
well as by heavy propagators. Which decay mode 
dominates depends strongly on the parameters at the 
GUT scale. We will discuss in detail this topic in the 
next section. 

As a third type of B3 mSUGRA scenarios we want to 
mention t\ LSP scenarios with a dominant Ag^^. cou- 
pling. Here, the dominant B3 operator couples directly 
to the fi LSP and allows for a 2-body decay of the t\ 
into two jets, 



A',, 



Uj,dk . 



(III.44) 



Agjj, can not generate A333 via the RGEs, because Xijk 
has to be anti-symmetric in the indices i, j. A3„„ with 
71 7^ 3 will be generated by the muon (n = 2) or elec- 
tron (n = 1) Higgs Yukawa coupling, c/. Eq. (jll.32[) . 
But since these Yukawa couplings are so small, the 
decay via A3„„ is too small to be seen. 

For j = 3; the up- type quark in Eq. (jIII.44p is a top 
quark and hence the decay Eq. (jIII.44p is kinematically 
forbidden for m^^ < mt- The fi LSP than decays in 
a 3-body decay mode via a virtual top quark into a 
W boson and two jets, where at least one jet is a & jet. 



W-bdk 



(III.45) 
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FIG. 4: Feynman diagrams contributing to the 4-body decay t-^ ^ t jj, Ujdk of the fi LSP via \'2jk- In this example 
the fi decays via a virtual neutralino x? = li 2, 3, 4) into a tau r~, a muon an up-type quark Uj of generation j 
and a down-type anti-quark c?fe of generation k. 






FIG. 5: Feynman diagrams leading to the 2-body decays of 
the fi LSP via the generated coupling A233 • The fi decays 
either into a muon fj,^ and a neutrino or into a t~ and a 
neutrino. 



We present the squared matrix element and the partial 
width of this process in Appendix [B] which to our 
knowledge has not been given in the literature so far. 



B. Dependence of ri Decays on mSUGRA 
Parameters 



In this section, we investigate the conditions at the 
GUT scale that lead to 2-body decays of the fi LSP. 
We assume a non- vanishing X'^^f, coupling at the GUT 
scale. This can easily be generalized to A'^^^. . We point 
out that the branching ratios of the fi LSP do not 
depend on the magnitude of A-^j,, since they cancel in 
the ratio. The following discussion is therefore also 
applicable to scenarios where the couplings are too 
small to produce a significant number of single slcpton 
events at the LHC but where the fi LSP is produced 
in cascade decays of pair produced SUSY particles. 

For the numerical implementation we use Softsusy 
2.0.10 [1^ to calculate the mass spectrum at the 
SUSY scale, Eq. (|II.37|) . In addition, we use our own 
program to calculat e A^^^ and Ai33 at the SUSY scale 
as described in Sect. Ill El We than pipe the mass spec- 
trum and the couplings through Isawig 1 . 200, which 
is linked to Isajet 7.75 [931. Isajet calculates the 
2-body partial width of the SUSY particles and pro- 




FIG. 6; 2-body decay branching ratio as a function of tan (3 
for different dominating \'2jk couplings at the GUT scale. 
The quark mixing is in the down sector and the mSUGRA 
parameters are Mo = GeV, M1/2 = 500 GeV, Ao = 
600 GeV, sgn(^t) = +1. 



duces an output for Herwig [981, l99|, |100|. We use a 
special version of Herwig 6 . 510 which also calculates 
the 4-body decays of the fi LSP [l06j . As an output, 
we consider the total 2-body decay branching ratio of 
the fi LSP, BR2. It is defined as 



BR2 = 



1 



1 + r4/r2 ' 



(III.46) 



where r2 and r4 denote the sums of the partial widths 
for the 2- and 4-body decays, respectively. 

We first show in Fig. [6] (Fig. [7|) the tan/3 dependence 
of the 2-body decay branching ratio. We give values 
for different non-vanishing couplings Ajj-j, at the GUT 
scale and we assume quark mixing in the down (up) 
sector. 

Nearly all fi LSPs will decay via a 2-body decay for 
large values of tan/?, i.e. tan/3 > 30, and down-type 
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tan/3 Ao [GeV] 



FIG. 7: 2-body decay branching ratio as a function of tan P 
for different dominating \'2jk couplings at the GUT scale. 
The quark mixing is in the up sector and the mSUGRA 
parameters are Mo = GeV, M1/2 = 500 GeV, A. 







600 GeV, sgn(/x) = +1. Couplings \'2jk for which the 2- 
body decay branching ratio nearly vanishes are not shown. 



mixing. In the case of up-type mixing this is also true 
for A211, A222 a-iid '^233- This behavior can be easily 
explained with the help of Eq. (|III.46p . The partial 
widths r4 can be approximated by (lo| 



T2 oc a: 



233 ' 



Ta oc A: 



X f 



(III.47) 
(III.48) 



TO^ denotes the mass of the relevant gaugino and m j 
denotes the mass of the virtual sfermion which couples 
directly to L2QjDk, cf. Fig.gl 

As we argued in Sect. Ill El the generated couphng A233 
scales roughly with tan^/3, cf. Eq. (jll.39p . Therefore, 
r2 scales with tan*/3. At the same time, X211 is hardly 
affected by tan /?. This is the main effect that enhances 
BR2 for large tan /3. 

Furthermore, increasing tan f3 increases the contribu- 
tion from the tau Yukawa couplings to the various 
RGEs. This is encoded in the function X^., Eq. pi.l3|) 
which is proportional to (l-|-tan^ /3). As can be seen in 
Eq. pi.lSp . increasing tan/3 and Xr reduces the mass 
of the right- and left-handed stau and therefore, with 
Eq. (|II.14p . the mass of the fi LSP, mf^. Furthermore, 
the off-diagonal matrix elements of the stau mass ma- 
trix Eq. (jlLlip also increase with tan (3. This leads to 
a stronger mixing between the right- and left-handed 
stau and lowers the mass of the fi, cf. Eq. pi.l4p . 

Note that r4/r2 is proportional to m| . According to 



FIG. 8: 2-body decay branching ratio as a function of Ao 
for non vanishing A211 at the GUT scale and different tan j3. 
We assume down-type quark mixing. The other mSUGRA 
parameters are Mo = GeV, M1/2 = 500 GeV, sgn(/i) = 
+f. The solid red curve corresponds to tan/3 = 7. 



Eq. (jIII.46p . the 2-body decay branching ratio there- 
fore strongly increases for decreasing rrif-^ . 

We observe in Fig. [6] also a large hierarchy between 
the different couplings X'2jk- example, a domi- 

nant A233 coupling leads to BR2 ~ 100% for any value 
of tan/3, whereas for X211 this is only the case for 
tan/3 > 25. This hierarchy reflects the hierarchy of 
the down quark Yukawa matrix elements, Eq. (jll.38p . 
which enter as the dominant term in the RGE of A233, 
Eq. (ITL321) . 

For up- type quark mixing. Fig. [71 and j ^ k the down- 
quark Yukawa matrix elements and therefore BR2 are 
nearly vanishing. 

We investigate the dependence of BR2 on Aq in Fig. [H 
for a dominant coupling A211 and down-type mixing. 
We see a minimum at Aq w 250 GeV. Here, BR2 is re- 
duced by up to 70% compared to Aq — ±1 TeV. The 
minimum and the position of the minimum is domi- 
nated by the following two effects. 

The right-handed stau couples to a left-handed stau 
(tau sneutrino) and a neutral Higgs (charged Higgs) 
via a trilinear scalar interaction {hE)33 [lOj . The cou- 
pling (h£;)33 has dimension one and in mSUGRA mod- 
els it is equal to Aq x (Y£;)33 at the GUT scale. The 
RGE of the right-handed scalar tau mass, m^^, de- 
pends in the following way on (h£;)§3 [lo| : 



^ 

dt 



(III.49) 



This term decreases m?^ when we go from the GUT 
scale to the SUSY scale pi.37|) due to the plus sign. 
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FIG. 9: 2-body decay branching ratio as a function of Adi/2 
for non vanishing A211 at the GUT scale and different tan /3. 
We assume quark mixing in the down sector. The other 
mSUGRA parameters are Mo = GeV, Ao = 600 GeV, 
sgn(/^) = +1. The solid red curve corresponds to tan/3 = 7. 



FIG. 10: 2-body decay branching ratio as a function of 
Mo for non vanishing A211 at the GUT scale and different 
tan/3. We assume quark mixing in the down sector. The 
other mSUGRA parameters are M1/2 = 1400 GeV, Ao = 
600 GeV, sgn(/i) = +1. 



The (negative) contribution of this term to m?^ is 
proportional to the integral of (Iie)'^-^ from tmin — 
ln(/V/GUT) to tmax = In(Mz). For the mSUGRA 
parameters given in Fig. [U Mq — GeV, M1/2 — 
500 GeV, sgn(/i) = +1, the integral of {hE)^^ is min- 
imal at ^0 ~ 180 GeV and, therefore, mf„ is maxi- 



this also leads to a maximum of 



rrv^^ and hence to a minimum of BR2 . 



mal. For 

r4/r2 

But the lightest stau is an admixture of the right- 
and left-handed stau. The off-diagonal mass matrix 
elements Blr, Eq. (jll.lip . depend also on the value 
of (hE)33 at the SUSY scale, Eq. ((IL37)) . through 
Ar = (h£)33/(YB)33. For ^0 = 180 GeV we find 
At Ki —110 GeV. A negative value of Ar enhances the 
effect of L-R-mixing which decreases mf ^ . Therefore, 
the maximum of as a function of Aq is shifted 
to Aq ~ 250 GeV compared to m^^. Note however 
that the At dependence of stau L~R-mixing is sub- 
dominant around the minimum because of /itan/3 ^ 

Ar. 

Next, we study the dependence of BR2 on the uni- 
versal gaugino mass Mi/2- We show this behavior in 
Fig. [HI again for a dominant X211 and down-type mix- 
ing. The 2-body decay branching ratios approach a 
constant value for increasing Mj^/2- Both, the squared 
mass of the gauginos, cf. Eq. (|II.15p . and the squared 
masses of the sfermions, cf. Eq. pi.lOp . depend lin- 
early on 

Aff/2- Therefore, 



M- 



lim r4/r2 oc 



? 4 



constant . 



(III.50) 



The dependence of BR2 on M1/2 for A/1/2 < 1 TeV is 
more involved, because the ratio r4/r2 depends also 
on the other mSUGRA parameters, mainly through 
the running sfermion masses, cf. Eq. pi.lOp . For ex- 
ample, we observe in Fig. ^ that the slope of BR2 
for Ml/2 ^ 1 TeV strongly depends on tan/3. For 
tan/3 =10, the slope is small and positive whereas for 
tan/3 > 13 the slope is negative. The magnitude of 
the slope also increases when we consider larger val- 
ues of tan/3. This behavior is again related to the 
tau Yukawa coupling (Y£;)33 and its effects on the fi 
mass described by the function Xr, Eq. pi.lSp . For 
large values of M1/2, the influence of Xr on the fi 
mass nearly vanishes. But as we go to smaller values 
of Mi/2 the (negative) contributions due to (Y£;)33 
become more and more important. For example, for 
tan/3 = 22 and Afi/2 = 1 TeV (M1/2 = 400 GeV) the 
Xr term reduces the mass of the right-handed stau 
by 3% (10%) compared to vanishing (Y_e)33- This re- 
duction of TTifj will also reduce r4/r2 resulting in an 
increase of BR2. This effect is more pronounced for 
large tan/3 because Xr is proportional to (1 -|-tan^ /3). 
If we neglect the effect of {Ye)33, the BR2 curves in 
Fig. [2]all get a small positive slope. 

Finally, we show in Fig. (TU] the dependence of BR2 
on the universal soft breaking scalar mass Mq. Here, 
we have chosen a rather large value of M1/2, A/1/2 = 
1400 GeV, because otherwise a fi LSP would exist 
only in a small interval of A/g. 

The behavior of BR2 can easily be understood. In- 
creasing Mq increases the mass of the sfermions, 
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and second to discuss the special cases Aj^j, ^ 0|gut 



and X'^jj, ^ 0|gut in more detail (Sects. IIVBI and 
HVCll . 



A. General Signatures 



FIG. 11: Example Feynman graph for single slepton pro- 
duction in fi LSP scenarios where the slepton decay 
proceeds via the generated A233 coupling (2-body decay 
mode). 



Ml 




FIG. 12: Example Feynman graph for single slepton pro- 
duction in fi LSP scenarios where the slepton decay pro- 
ceeds via the dominant \'2j^. coupling (4-body decay mode). 



Eq. pi.lO|) . but not the mass of the gauginos. There- 
fore, the nominator of r4/r2 oc TO|^/(m|m^) is a poly- 
nomial of order 0{Mq ), whereas the denominator is 
only a polynomial of order C'(Mq). Therefore, the 2- 
body decay branching ratios fall off for increasing A/q 
as shown in Fig. [TOl The lines in the Figure termi- 
nate at values of A/q above which the fi is no longer 
the LSP. 



IV. RESONANT SINGLE SLEPTON 
PRODUCTION IN n LSP SCENARIOS 

We now apply the previous discussion to resonant 
single slepton production in B3 mSUGRA scenarios 
with a fi LSP. Charged sleptons ^Li and sneutrinos Vi 
can be produced singly on resonance at the LHC via 
qk qj annihilation processes. The production cross sec- 
tion is proportional to |A-^j.p and therefore large slep- 
ton production rates are expected in scenarios with a 
dominant A^^ j. coupling. The RGE generation of A^aa is 
important for the subsequent slepton decay in fi LSP 
scenarios. As discussed in the previous section, a non- 
vanishing Ai33 introduces new 2-body decay channels 
for the fi LSP. The interplay of these 2-body decays 
and the 4-body decays via A'^^^ determines the final 
state signatures. In Figs. [TT] and [T^l example Feyn- 
man graphs for single slepton production and the sub- 
sequent decay in fi LSP scenarios are shown. 

It is the aim of this section to first give a general 
overview of the possible final states for these reactions 



In the last section, the ratio of 2- to 4-body fi LSP 
decay rates and its dependence on various SUSY pa- 
rameters has been studied. Now, we focus on single 
slepton production in fi LSP scenarios and are inter- 
ested in the general decay patterns, independent of the 
precise SUSY parameters. We first give an overview 
over all possible final states and signatures which could 
be used as the starting point for an experimental anal- 
ysis. 

A (left-handed) charged slepton or sneutrino can be 



produced directly via 
modes: 



X'^j). and has several decay 



"■Li 



dj dk 



Am 7 



0+ 

n 7 



(IV.51) 



(IV.52) 



Both can decay via the B3 coupling, which is the in- 
verse production process. It is however suppressed by 
l^ijfeP- If Kjk ^ O{10-^), it contributes typically at 
the percent level. The dominant decay channels are 
2-body decays into a lepton-gaugino pair. Further 3- 
and more-body decays are expected to be negligible, 
due to phase space suppression. 

In case of j = 3, the hadronic production of a charged 
slepton cannot proceed via two quarks as given in 
Eq. (jlV.Sip . due to the vanishing top-quark parton 
density inside a proton. Instead, the slepton can for 
example be produced via a gdk initiated Compton pro- 
cess in association with a single top quark. Further- 
more, the decay into tdk may be kinematically for- 
bidden. In this case, the slepton decays via a vir- 
tual top. The corresponding decay width is given in 
Appendix |B] Sneutrino production for j = 3 is pos- 
sible, Eq. (jIV.52p . but due to the low bottom-quark 
density small cross sections are expected. We do not 
consider 7 = 3 any further here and refer the reader 
to [53, [63, [zOj [zII for a detailed investigation of this 
topic. 

For the following discussion, we assume that the pro- 
duced slepton predominantly decays into a lepton and 
the lightest neutralino. This assumption is motivated 
by the fact that we consider t\ LSP scenarios. In these 
scenarios, sleptons are light compared to gauginos and 
decays into heavier neutralinos or charginos will be 
kinematically excluded or strongly suppressed. See 
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also the computed branching ratios in explicit SUSY 
models in 

The produced Xi is not the lightest SUSY particle and 
will decay further into the fi LSP, 

Xi^r^ft- (IV.53) 

Since the neutralino is a Majorana fermion, both 
charge conjugated decays are possible. In most fi LSP 
scenarios this is the only possible decay mode of the 
neutralino. However, in some scenarios, the right- 
handed sleptons fiR and are lighter than the Xi 
and the additional channels Xi — *■ ^r^^ are open (for 
£ = fj,,e). The £r subsequently decays into the fi LSP, 
a r, and a lepton via a virtual neutralino 

These decay chains have smaller BRs than the decays 
in Eq. pV.53|) . However, they lead to an additional 
lepton pair in the final state and could be, therefore, 
of special interest for experimental analyses. 

B. A'ajfe / 0|gut, A233 < y2jk 

Let us now study more detailed the final state sig- 
natures in a scenario with Aj^j, 7^ 0|gut and a gener- 
ated A233 coupling which is small but non-zero at lower 
scales. In these scenarios, resonant single (11 produc- 
tion and resonant single production at hadron col- 
liders is possible, 

Uj dk ^ fll ^ Uj dk/pL" Xi, 

(IV.55) 

dj dk i>f_,^ dj dk/Vf^ Xi- 

As explained above, a small fraction of the sleptons 
decay via the inverse production process. Predomi- 
nantly they decay into a lepton and the lightest neu- 
tralino, Xi- The decays involving heavier neutralinos 
or charginos are typically not accessible. 

The difference between ji l and production concerns 
the flavor of the initial quarks involved (which is re- 
lated to different parton density functions and is thus 
important for the hadronic cross sections), and the 
nature of the lepton resulting from the slepton decay. 
In both processes a neutralino is produced in the pre- 
dominant decay, which in turn decays into the fi LSP, 
as given in Eq. pV.53p and Eq. (|IV.54p . Finally, the 
fi decays either via the dominant X2ji. coupling (4- 
body decay) or via the generated A233 coupling (2- 
body decay). For the 4-body decays, only the decays 
via virtual neutralinos have to be considered. Decay 
modes via virtual charginos are suppressed due to the 
larger mass and their weaker couplings to the predom- 
inantly right-handed fi LSP. The complete cascade 
decay chains are listed in Tab. IIIII 





Uj dk ^ ti^ > Uj dk/fi Xi 

or 

dj dk — > dj dk/vp.xR 


Y? r+ rr 
[x^i^r^r-t+r] 


Y? ^ T- f,+ 
Al ' ' 1 




— > T^/i" Uj dk 
ff — > fl'^Ujdk 

f[' T~v^ dj dk 
T^v^ dj dk 


f + T^fl'^ Uj dk 

f + -r+^" Uj dk 
T^iy^ dj dk 
T^y^i dj dk 


A233 




ft r+u^ 



TABLE III: Slepton decay chains with all possible final 
states for single jl^ and single production via X'2jkj 
respectively. The charge conjugated processes are not 
shown explicitly. Slepton decays into heavier neutralinos 
or charginos are neglected. The Xi decays predominantly 
into a fi LSP and a r. In some scenarios, decays as in 
Eq. HIV.54P are possible, they are cited in brackets. Owing 
to the Majorana type nature of the neutralino two charge 
conjugated decays of the neutrahno are possible (second 
and third column) . In the first column the B3 coupfing in- 
volved in the subsequent 4- or 2-body fi decays are given. 



A classification of all possible final state signatures is 
given in Tab. IIVI for //^ and for production. For 
completeness, we include here the direct B3 decays via 
Ajjj,, which usually contribute at the percent level for 
couplings at the order of 0(10"^). Neutrinos do not 
give a signal in a detector and are denoted as missing 
transverse energy, ^t- Final state quarks are treated 
as indistinguishable jets, j. 

The 4-body decays via Xi^j/. and the 2-body decays via 
the inverse production process lead to two jets in the 
final state. In contrast, the 2-body decays via A233 are 
purely leptonic. Many cascade decay chains provide 
missing transverse energy. Furthermore, since we are 
considering fi LSP scenarios, there is always at least 
one T among the final state particles. The experimen- 
tally most promising signatures are most likely those 
involving a large number of muons, for example like- 
sign dimuons and three or four final state muons. If 
the Xi decays only into rir, there are two signatures 
including like-sign dimuons for //^ production. For 
i/^ production, muons can be produced singly only. 
But if the decays Eq. pV.54[) are open, both slepton 
production processes allow for dimuon and trimuon 
production. In case of p,L production, even four final 
state muons are possible. Additionally, depending on 
how easily taus will be identified, an analysis could be 
based on like-sign /ir-pairs. 

The final state signatures depend sensitively on which 
particle is the LSP. Compared to slepton production in 
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the xi LSP scenarios [H, [H, [H, [H, [H, H [U, [H, [H , 

there are three main differences here. First, for a 
fi LSP we have always one or two taus in the fi- 
nal state, which in x1 LSP scenarios is only possi- 
ble for smuon production if heavier neutralinos are in- 
volved in the decay chain. These heavy neutralinos 
then decay into the lightest ncutralino and possibly 
taus. Second, the generation of a A coupling can be 
neglected in Xi LSP scenarios. As argued above, A 
only allows for additional 3-body decays which arc 
thus not phase-space enhanced compared to the 3- 
body decays via the dominant A' coupling. As a con- 
sequence, purely leptonic final state signatures are ab- 
sent in Xi LSP scenarios. Third, due to the modified 
spectra in Xi LSP scenarios, also j/^ production can 
provide like-sign diniuon events. In this case, t'^ can 
often decay into a fi and a chargino. Like-sign dimuons 
arise cither if the chargino directly decays via A' into 
a n and two quarks, or if the chargino first decays into 
the Xi LSP and then the Xi LSP decays via A' into a 
fi and two quarks. 

This discussion can easily be translated to scenarios 
with A'j^^j, 7^ by replacing the muons by electrons 
(and vice versa). Since there is typically no differ- 
ence in mass between sleptons of the first and second 
generation, respectively, the kinematics arc the same. 
Note however that the bounds on the B3 couplings are 
stronger for Aj^j, than for Aj^;. for example due to the 
non-observation of neutrinoless double beta decays. 



C. X'sjk / 0|gut 

Some additional remarks are in order for a dominant 
^'sjk ^3 coupling. These couplings allow for resonant 
single i>r production and, owing to the L-R-mixing in 
the stau-sector, also both resonant fi and T2 produc- 
tion (jV 3). 

For fi production, we refer to the discussion of LSP 
decay modes in Sect. IIII Al Here the LSP couples di- 
rectly to the B3 operator and the inverse production 
process dominates the decay rate, 

Uj dk f[' Uj dk ■ (IV. 56) 

This decay is kinematically accessible if j 7^ 3. For 
J = 3 the stau decays via a virtual top-quark, cf. 
Eq. pil.45p . for m^^ < mj. Note that j = 3 requires 
associated production, e.g. g dk — ^ f due to the ab- 
sence of top quarks inside the proton [sO, IH, [zO, [Zll • 

For f2 and production, there are the following 2- 
body decay modes: 

{Uj dk, 
T-x\ , (IV.57) 





jlj^ production 




r+ T" fi' fi+ jj 
r+T' 11' jj 
[r+r' fi' fi- fi+ fi+ jj] 
[r+r- fi~ fi- fi+ jj] 
[r+r' ft' fi+ e+e" jj] 
[r+r' fi' e+e- jj] 


A233 


, .-^ .,-f IT' 

t"*" t~ //~ 

[r± fi-fi-fi+fi+ ] 
[r+r- fi-fi-fi+ Pt ] 
[r+ ft- fi+ e+e- ] 
[r+r- ft- e+e- ] 


inv. prod. 


jj 




i'f_i production 


\' 


r+r- fi+ jj 
r+r- jj 
[r+r- fi- fi+ fi+ Pt jj] 
[r+r- fi- fi+ Pt jj] 
[r+r- ft+ e+e- jj] 
[r+r- e+e- jj] 


A233 


r ft^ Pt 

r+ r- 

[r+ fi-fi+fi+ ] 
[r+r- fi- fi+ ] 
[r+ fi+ e+e- ] 
[r+r- e+e- Pt ] 


inv. prod. 


jj 



TABLE IV; Summary of all possible final states for single 
slepton production via A2jj.. Decays involving the domi- 
nant A2jfc coupling and involving the generated A233 cou- 
pling are listed separately, cf. Tab. IIIII If kinematically 
allowed, the Xi m^-y a-lso decay into a light-flavor lepton- 
slepton pair which gives rise to an additional fi+ ft- or e+e~ 
pair in the final state. The corresponding signatures are 
given in brackets. The decay via the inverse production 
process is also listed. 



{dj dk, 
v.,x\ ■ (IV.58) 

The inverse production process contributes and leads 
to a jj final state. The decay into a lepton and a neu- 
tralino often dominates for small tan/3 (tan/3 < 10). 
The neutralino decays further into the t\ LSP which 
directly decays into two quarks: 

X;^r±f+, f^^u.dk, (IV.59) 

where we have included the two charge conjugated de- 
cays of the neutralino. The final states of these de- 
cay modes are T~r^Jj, and there is the possibility of 
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like-sign tau events. If the Xi decay (|IV.54[) is kinc- 
matically allowed, we can have an additional pair of 
electrons or muons in the final state. 

The singly produced slcpton can also decay into the 
fi LSP and a SM particle, h^, or W, respectively 
(final states: h'^ /Z^ /W jj). This decay mode is special 
for singly produced sleptons of the third generation 
because they are L-R mixed eigenstates. It can be the 
dominant decay mode of the f2 and Vr, depending on 
the parameters. 

The branching ratios for all B3 conserving T2 and i/r 
2-body decay modes are given in Tab. lXIl in Appendix 
iBl for the SUSY parameter sets A and B. 



V. SINGLE SMUON PRODUCTION: AN 
EXPLICIT NUMERICAL EXAMPLE 

In this section, wc present explicit calculations of 
promising signal rates for resonant slcpton production 
at the LHC in the B3 mSUGRA model with a fi LSP, 
focussing on parameter sets A and B, cf. Eq. (|II.16|) . 
First, we consider in Sect. IV Al (exclusive) like-sign 
dimuon events, i.e. events with exactly two muons of 
the same charge in the final state. An analysis of SM 
and SUSY backgrounds for the like-sign dimuon sig- 
nature is given in Sect. IV 51 Second, in Sect. IV Cl we 
present event rates for single smuon production lead- 
ing to three or four muons in the final states, which 
are kincmatically accessible within sets A and B. 



A. Like-Sign Dimuon Events 

Following Refs. [13, H^, we first concentrate on events 
with exclusive like-sign dimuons. Here events with 
more than two muons are rejected. In this sense, in 
fi LSP scenarios, only single smuon production leads 
to exclusive like-sign dimuon pairs, cf. Tab. IIVI It 
has been shown in Refs. [H, [5^ that this selection 
criterion enhances the signal to background ratio con- 
siderably. In Refs. [5^. Issj it was shown that using 
a set of cuts, the SM background rate at the LHC, 
TbIsMj can be reduced to 



S|SM 



4.9±1.6evcnts/10fb \ 



(V.60) 



At the same time the cut efficiency, i.e. the number 
of signal events which pass the cuts, lies roughly be- 
tween 20% and 30%. Note that Refs. [H,!!^ assume a 
Xi LSP. As we will argue in Sect. IV 51 similar cuts are 
also applicable in fi LSP scenarios. For the numbers 
presented in this section, however, no cuts are applied 
and full cross sections and event rates are given. 

The total cross section for like-sign dimuon events is 
given by the resonant /ij or jij^ production cross sec- 
tion multiplied by the respective branching ratios lead- 
ing to like-sign dimuon final states. Both decays via 



the dominant X'2ji; coupling and a generated A233 cou- 
pling contribute. For a negatively charged smuon they 



A' 



^1 






A' 






^+ 


T 




A 












A' 




Uj 
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plus the analogous decay chains where the neutralino 
decays first into an e^-e^ pair, cf. Eq. (|IV.54p . The 
couplings depicted on the arrows indicate the cm- 
ployed B3 coupling. The decay chain for a positively 
charged smuon can be obtained by charge conjugation. 
However, one should keep in mind that the production 
cross sections for fi^ and /i£ differ at pp colliders, since 
charge conjugated quarks (and corresponding parton 
densities) arc involved. 

The cross sections for the exclusive like-sign dimuon 
final states are presented in Tab. |V| for Set A and in 
Tab. IVII for Set B. The smuon production cross sec- 
tions, (Tprod.(/iL) (see also Tabs. HXl and IXl). include 
NLO QCD and SUSY-QCD corrections [gj, see Ap- 
pendix]^ For the numerical analysis, we only consider 
couplings \'2jk that involve partons of the first gener- 
ation leading to large production cross sections at the 
LHC. 

As already discussed, the fi LSP can cither decay 
via A' (4-body decay) or via A (2-body decay). A 
list of the respective branching ratios is given in Ap- 
pendix A, Tabs. IXIII and IXIIII for sets A and B and 
for several Aj^j. couplings. Here we show the result- 
ing cross section times branching ratio, CTprod. x BR^/ 
and CTprod. X BRa, for like-sign dimuon events involv- 
ing fi decays via A' and A, respectively, as described 
in Eq. (fOTjl . 

The total number of exclusive like-sign dimuon events 
is given by the integrated luminosity multiplied by the 
total cross section. In Set A with up-type (down-type) 
quark mixing, we obtain per lGfb~^ 



A^(Ai"Ai" +Ai+^+)/10fb" 

prod.(AL) + CTprod. (Ml) 



'325 (330) 
110 (115) 
195 (210) 

,110 (115) 



BRa' + BRa 

^^211 



X 10 



/lOfb"^ for 



-^221 
^^212 
-^213 



0.002|gut. 



(V.62) 



Note that for up-type mixing, some larger couplings 
may be considered. From the neutrino mass bounds, 
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up-type mixing down-type mixing 
IT J X RR X / rr J X RR ^ IT J X RR , , n j x RR % 

Oprod. ^ J^J^A' t-'prod. ^ J^J^A t^prod. ^ J-*J^A' t-'prod. ^ J->J^A 


A^n = 2 X 10-^ GUT + + 

/i 108 


11.1 0.71 9.81 2.09 

n A IOC 1^0 o/^/^ 

19. 4 1.25 17.2 0.66 


, 3 fi- 42.0 

A221 = 2 X 10 Igut + + 

fl fl 10. Z 


7.84 - 4.51 3.88 

Q riQ 1 7/1 1 P^n 
o.Uo — 1.(4 I.OU 


, M~ 18-6 

A^12 = 2 X 10-3 GUT + + 

/I /i^ 86.0 


3.46 - 1.99 1.71 
16.1 - 9.23 7.94 


3 ^"M" 8.80 
A2i3=2xl0 Igut ^+ ^+ 


1.67 - 1.32 0.40 
9.43 - 7.43 2.24 



TABLE V: Cross sections for exclusive like-sign dimuon (fi^ fi^ or final states at the LHC within Set A. In the left 

column, we present the single-smuon production cross sections, (Tpiod.(/iJ), see also Tabs. HXl and 1x1 In the right column, 
we have folded in the relevant decay branching ratios, in order to obtain like-sign dimuons. All cross sections are given in 
fb. Where they exist, we have assumed always a cascade of 2-body decays. We consider in turn quark mixing in the up- 
and down-sector, when determining the dominant fi decay mode. The fi LSP can either decay via A' (4-body decay) or 
via A (2-body decay), cf. Tab. IIIII which leads to different like-sign dimuon cross sections, cTprod. x BR;^/ and aprod. x BRa, 
respectively. The A2jfc couplings are in accordance with neutrino mass bounds [13, [121 ■ In case of up-type mixing, larger 
values of A2jfe for the four considered couplings are allowed by the neutrino mass bounds. The cross sections scale with 
|A'p and the corresponding rescaling can easily be performed. 



^ CTprod. (Al) [fb] 


up-type mixing down-type mixing 

O-prod. X BRx' O-prod. X BRa O-prod. X BRx' CTprod. X BRa 


„, M~ 476 

A^. = 1 X 10-lGUT ggg 


1.04 101 0.21 102 
1.93 188 0.39 189 


, 2i M /^^ 309 

A221 = 1 X 10 GUT + + 

fj, fi 105 


62.8 - - 66.2 
21.4 - - 22.5 


/ 2i M /^^ 123 

A212 = 1 X 10 GUT + + 

fj, fi 681 


25.1 - - 26.3 

139 - - 146 


, 2i M /^^ 54.6 
A213 = 1 X 10- Igut ^+ ^+ 


11.2 - 0.02 11.7 
75.6 - 0.16 79.4 



TABLE VI: Same as Tab.FVlbut for single slepton production within Set B. The neutrino mass bounds are less restrictive 
in the case of Set B and Xi^ji^ — 0.01|gut are considered for both up- and down-type quark mixing. All cross sections are 
given in fb. 



also \2ii 221.212.213 = 0.01|gut (and even larger) are 
allowed. The cross sections arc proportional to |A'p 
and thus a five times larger coupling implies cross sec- 
tions and event numbers multiplied by a factor of 25 
compared to those of Tab. |Vl 



For Set B, A2J7,. = 0.01|gut is allowed for both up- and 
down-type mixing. The numbers of like-sign dimuon 



events are, 

N{^i-^i- + M+^+)/10flo-^ 

'2920(2920) 
840 (890) 
1640(1720) 
t870 (910) 



^211 



/lOfb-^ for 



A' 



221 



^212 



^213 



= 0.01|gut, 



(V.63) 



for up-type (down-type) quark mixing, respectively. 

As can be seen in Eqs. (|V.62p . (|V.63|) . for each non- 
zero A' coupling the total event numbers for up- and 
down-mixing are of the same order. But as Tabs. fVl 
and IVII show, the parts contributing to the event rate 
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can be quite different. In case of up- type mixing and 
j k, the 4-body decays via A' dominate and the 
contributions of the 2-body decay are ncgUgible [since 
the size of the necessary A coupUng is proportional to 
(Yu)^^]. In contrast, for down- type mixing all four 
considered couplings can generate a relatively large 
A233, cf. Fig. [21 and the 2-body decay modes con- 
tribute considerably. In Set B. where tan/3 is large 
and where thus the fraction of 2-body decays is es- 
pecially high (see discussion of Fig. [5]), reliable event 
numbers are only obtained if the generation of A233 
is included in the theoretical framework. Moreover, 
a measurement of the ratio of 2-body to 4-body fi 
decays can reveal information about where the quark 
mixing takes place. 

For j = k, the generation of a A coupling is also possi- 
ble in case of up-type mixing. In Set A, the generated 
A233 is not large enough to allow for large 2-body de- 
cay rates. However in Set B, due to the large tan/3 
value, the 2-body decays dominate over the 4-body 
decays. Thus, the different fi decay modes contain 
also information about tan /3. 

We present in Tabs. fVl and IVII also the total hadronic 
cross sections for single smuon production, crprod.(Mi)- 
Within one parameter set, the cross sections vary 
strongly for different Aj^j,. This is of course related to 
corresponding required parton density functions. The 
largest cross section is obtained for A'j^i 7^ 0, i.e. for 
the processes ud ^ p.^ and ud ^ jX^. Smaller cross 
sections are obtained for A212 7^ (involving an up 
quark and a strange quark) and the smallest cross sec- 
tion for A221 7^ (charm quark and down quark) and 
•^213 7^ (up quark together with bottom quark). 

Since the LHC is a pp collider, there is an asymmetry 
between the /ij and jT^ production cross sections. If 
experimentally a distinction between fJ.^ fJ-'^ and fJ.~ fJ.~ 
event rates is found, the ratio can be used to constrain 
the indices of the non-zero Aj^^. coupling. For exam- 
ple, a non-vanishing coupling X'211 leads to a ratio of 
N{iJ,^fj,~^) : N{iJi~ iJL~) ~ 2 : 1 in sets A and B, whereas 
for non- vanishing A221 the ratio is 1 : 2.5 in Set A and 
1 : 3 in Set B. The highest event rates are obtained 
for processes that involve the valence quarks u and d. 
The charge conjugated processes, involving u or d, are 
suppressed in comparison. Thus, a larger fraction of 
fi'^ fj-'^ events goes along with j — 1 (where the pro- 
duction process is u d^ p.^ ) and a larger fraction of 
IJ,~li~ events is related to fc = 1 and j 7^ 1 (production 
process ujd ^ p.^)- 

B. Discussion of Background and Cuts for 
Like-Sign Dimuon Final States 

In this section, we discuss the background for like-sign 
dimuon events from the SM and from SUSY parti- 
cle pair production via gauge interactions. We follow 
Refs. [52I, [531 closely. There, single smuon production 



via A211 was investigated assuming a LSP. A de- 
tailed signal over background analysis was performed 
based on like-sign dimuon events. We argue that a 
similar or even the same set of cuts might be used to 
suppress the background in our case and we compare 
background and signal rates to determine the discov- 
ery potential of our analysis. 

The main SM background sources are ti production, 
bb production, single top production, and gauge boson 
pair production, i.e. WW, WZ and ZZ production. 
In Refs. (5^. [ssj. the dominant signature from single 
smuon production including like-sign dimuon events is 

Pl^^l-x1^^^-{^i-ud). (v.64) 

The two muons of the signal (|V.64p are isolated be- 
cause they stem from different decays of SUSY par- 
ticles. In addition, the muons carry large momenta 
since they originate from the decay of (heavy) SUSY 
particles. The following cuts were proposed to improve 
the signal over SM background ratio at the LHC: 

• The muon rapidity |?7| < 2.0, thus requiring all 
the leptons in the central region of the detector, 

• a cut on the transverse momentum on each 
muon: — 40 GeV, 

• an isolation cut on each of the muons, 

• a cut on the transverse mass of each of the 
muons, 60 GcV < Mr < 85 GeV, 

• a veto on the presence of a muon with the oppo- 
site charge as the like-sign dimuons, 

• a cut on the missing transverse energy, ^ 
20 GeV . 

These cuts reduce the SM background to 4.9 ± 1.6 
events per 10 fb-^ at the LHC , cf. Eq. dEMI. Among 
the above cuts, the isolation and pt cut lead to the 
strongest suppression of the SM background. 

We now investigate the case of a fi LSP. If the 4-body 
decays (jIII.42[) of the fi LSP dominate, the leading 
signature of resonant single smuon production includ- 
ing like-sign dimuon events can be written as 

PI ^ /i-x! ^ ^l-T^f^ ^ ^l-T^iT^^l~ud). (v.65) 

As above, the muons originate from the decay of heavy 
particles (fi and Pl), are in general well isolated, 
and carry large momenta. Thus, for both signals 
Eq. (|V.64[) and Eq. (|V.65p . the same cuts should al- 
low to discriminate between the signal and the SM 
background. Furthermore, the additional pair of taus 
in Eq. (jV.65|) allows to require one or two (isolated!) 
taus. This might additionally improve the signal to 
background ratio. 
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If the fi LSP predominantly decays via 2-body decay 
modes, Eq. pil.43p . the situation is a bit different. 
The hke-sign dimuon signature is now 

Ml M^Xi /i"T+f" frr+ii^ri^r)- (V.66) 

We again have two isolated muons with large momenta 
and the same isolation and pt|/j cuts as before should 
be useful to suppress the SM background. But the 
neutrino of the fi decay leads to high missing trans- 
verse energy /J^ in the signal and an upper bound 
on is not appropriate anymore. Alternatively we 
propose a cut that requires a minimum missing en- 
ergy, e.g. f^T > 60 GeV. This would also reduce the 
SM background where the main source of are low- 
energetic neutrinos from W decays. Furthermore, we 
can again require an additional tau in the final state. 
Finally, one can exploit the fact that the 2-body de- 
cays lead to a pure leptonic final state and a jet veto 
can be applied. 

In Refs. [H, [11], the SUSY background on hke-sign 
dimuon events is suppressed by vetoing all events with 
more than two jets of prljet > 50 GeV. This cut will 
also work if the 4-body decay mode of the f\ LSP 
pil.42[) dominates. The 2-body decay modes lead to 
purely leptonic final states and even no high-p^ jet 
may be required. 

We conclude that for fi LSP scenarios, the background 
for like-sign dimuon events can be suppressed similarly 
as it has been proposed for x° LSP scenarios in [5^.[53|. 

We thus compare our signal, as given in Eq. (|V.62p and 
Eq. (|V.63p for sets A and B respectively, to the back- 
ground, assuming that cuts as discussed above reduce 
the SM background to less than 5 events per 10 fb~^, 
c/. Eq. (jV.60p . For the signal efficiency, we assume 
20%, i.e. 20% of signal events pass the cuts. We ne- 
glect systematic errors, at this stage of the analysis. 

For Set A a more than 5cr excess over the SM back- 
ground can be obtained for an integrated luminosity of 
10 fb~^ for all couplings given in Eq. (|V.62p . For Set 
B, a cut efficiency of 20% for the signal corresponds to 
an excess between 100 a and 300 a for the number of 
like-sign muon events over the SM background! There- 
fore, within Set B, couplings can be tested at the LHC 
down to \'2jk\GVT^ 0(10"^). But a detailed Monte- 
Carlo based signal over background analysis remains 
to be done. 



C. Final States with 3 and 4 Muons 

To round off our studies, we consider in this section 
final states with more than two muons. For exam- 
ple, for parameter sets A and B, the Xi cannot only 
decay into a fi-r pair but also into a jiR-y. or cr- 
e pair. These are kinematically accessible and have 
non-negligible branching ratios (Set A: 7.0%, Set B: 
2.2%; see Tab. |Xl|. As we have shown in Tab. |lVl 



these decays lead to three or even four muons of mixed 
signs in the final state. Each of the muons stems from 
the decay of a different SUSY particle. Especially the 
four-muon final state cannot be found at a high rate 
in x\ LSP scenarios and its observation could be a 
hint for a fi LSP. Therefore, we analyze the three- 
and four-muon final states in this section. All nec- 
essary branching ratios and production cross sections 
are given m the Appendix, see Tabs. HXllXlIIl 

The four-muon events may be classified into 
fi~ fi~ fj.^ , fi~^ ^'^ , and yU^/i+^+yU"'' signatures 

and we introduce the notations a{ h), (t( 1 — h), 

and (t(-|- + H — ), for the respective cross sections. The 
four-muon final states require a long decay chain and 
many different decays contribute at various stages. For 
smuon production, summing up all contributions, the 
cross sections can be written in the following compact 
form 

+) ='7prod.(/iZ) X BR(/i£ XifJ-^) 

+ +-) ='^/i( ^) X C^prod.(Mj)/CTprod.(/iZ) ' 

crr^i- - ++) = c^^l + ) + + +-)' 

(V.67) 

where Pf,{lfx) = BR(ff . . .) + BR{t+ 

. . .) denotes the probability of a negatively charged 
final state muon in a fi decay. The difference between 

ajx{ h) and o'^(-|--l--l — ) stems from the different 

partons and parton densities involved in the produc- 
tion cross sections. 

Smuon production can also lead to exactly three final 
state charged muons, n~ ^~ fX^ or jX^ . The cor- 
responding cross sections now involve the probability 
Pf J (0/i) for a fi decay without a final state muon, 

^) ='^prod.(Ai;) X BR(/i^ ^ X?M^) 

xBR(x?^m5,/^")x2P^,(Om), 

+ -) =0'a(- - +) X 0-prod.(Mi)/CTprod.(/iZ) ■ 

(V.68) 

There are 16 different decay chains of the /i^ leading to 
a iJL~ fX^ final state. The factor of 2 in Eq. (|V.68p is 
a consequence of summing over all these decay chains. 

The same final state signatures (exactly three muons) 
can be obtained via production. The decay chain 
is similar to that of a produced smuon. The missing 
muon from the slepton decay is here replaced by de- 
manding a muon in the final fi decay, 

Uui h) = [o-pi.od.(i>^) + crprod.(£'^)] 

X BR(:>^ X? i^f.) X BR(x? ^ Afl M") 
xPn(lM), 

<Jy{+ + -) =CT^(--+). 

(V.69) 
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Set A 








a(+ + +-) 


+) 


E^(-----) 


E^(+ + ...) 


-^211 


= 2 X 10"^ 


GUT 


9.38 (9.39) 


12.9 (13.0) 


5.32 (5.26) 


3.39 (3.35) 


1.93 (1.91) 


16.6 (16.6) 


21.7 (21.6) 


•^221 


= 2 X 10"^ 


GUT 


5.77 (5.77) 


3.84 (3.74) 


1.89 (1.77) 


0.53 (0.49) 


1.36 (1.27) 


9.02 (8.81) 


6.26 (6.00) 


-^212 


= 2 X 10"^ 


GUT 


4.02 (3.93) 


9.05 (9.24) 


3.39 (3.17) 


2.79 (2.61) 


0.60 (0.56) 


8.01 (7.66) 


15.2 (15.0) 


-^213 


= 2 X 10"^ 


GUT 


2.04 (2.02) 


5.14 (5.19) 


1.85 (1.80) 


1.57 (1.53) 


0.28 (0.27) 


4.17 (4.09) 


8.56 (8.52) 



TABLE VII: Cross sections for signals with three or four final state muons within parameter Set A, assuming down-type 
(up-type) quark mixing. Given are the cross sections as defined in Eqs. HV.67|) - (|V.70|I and the sums for two negatively or 
positively charged muons, ^ a{ . . . ) or ^ a{+ -I- . . . ), respectively. All cross sections are given in fb. 



Set B 


'^(--+) 


'^(+ + -) 


'^(--++) 


a(+ + +-) 


+) 


E^(-----) 


E'^(+ + ---) 


-^211 


= 1 X 10"^ 


GUT 


20.8 (20.8) 


29.1 (29.1) 


13.4 (13.4) 


8.73 (8.73) 


4.69 (4.69) 


38.9 (38.9) 


51.3 (51.3) 


•^221 


= 1 X 10"^ 


GUT 


11.9 (12.0) 


7.77 (7.59) 


4.08 (3.88) 


1.04 (0.98) 


3.05 (2.89) 


19.1 (18.7) 


12.9 (12.4) 


•^212 


= 1 X 10"^ 


GUT 


8.14 (7.98) 


19.5 (19.9) 


7.93 (7.53) 


6.72 (6.39) 


1.21 (1.15) 


17.3 (16.7) 


34.2 (33.8) 


^213 


= 1 X 10"^ 


GUT 


3.94 (3.85) 


10.4 (10.6) 


4.20 (4.00) 


3.66 (3.48) 


0.54 (0.51) 


8.68 (8.36) 


18.3 (18.1) 



TABLE VIII: Same as Tab. IVlJ but for single slepton production within Set B. All cross sections are given in fb. 



The total cross sections for (exactly) three final state 
muons are then given by 

a(TT±)=fTA(TT±) + ap(TT±). (V.70) 

Tabs. IVIII and IVIIII give an overview over the numer- 
ical results. The same A' couplings as in the previous 
Tabs.rVland lVII are considered. The generation of A233 
has been taken into account for the fi decays and the 
cross sections give total numbers, including both 4- 
and 2-body fi decays. 

We see that the sum of three- and four-muon events 
is in the same order of magnitude as the results for 
purely like-sign dimuons. For Set A, where BR(x? — > 
Mi? m) = 7%, the event numbers are even larger. In Set 
B, with BR(x]' flu fi) = 2%, the total contributions 
are smaller by a factor of about three. Depending 
on the experimental goals, these channels thus give 
important contributions and should be included in an 
analysis. On the other hand, these events also suggest 
to use three or four final state muons as a signal for 
slepton production since the background is expected 
to be very low. 

VI. CONCLUSION 

B3 interactions allow for LSP decays and thus reopen 
large regions in the SUSY parameter space, where the 
LSP is charged. We have investigated for the first time 
in detail the phenomenology of B3 mSUGRA models 
with a fi LSP. We have hereby assumed only one non- 
vanishing B3 coupling A-^j, at Mqut- 

An essential feature of the B3 mSUGRA signatures is 
the decay of the fi LSP. Given only one B3 coupling 



at A/gut, we would expect either a 4-body or 2-body 
decay of the fi LSP depending on whether it couples 
directly to the dominant B3 operator or not. However, 
in B3 mSUGRA models the RGBs are highly coupled 
and further couplings arc generated at the weak scale. 
These are of course suppressed relative to the domi- 
nant coupling but may lead to 2-body decays, which 
have larger phase space and do not involve heavy prop- 
agators. 

We have here numerically investigated the generation 
of Ai33 couplings via dominant A^^^ couplings. The 
generated couplings are typically smaller by at least 
two orders of magnitude; see Figs. [5] and O We have 
then performed a first detailed analysis of the param- 
eter dependence of the fi LSP decay modes. It turned 
out that in large regions of parameter space the 2-body 
decay dominates over the 4-body decay, see Figs.lMTUl 

In the second part of the paper, we applied our re- 
sults to resonant single slepton production at the LHC, 
which is possible in B3 scenarios with a non-zero 
A-^j, coupling. We first studied the general decay sig- 
natures. From the experimental point of view, the 
final states with two like-sign or even more charged 
leptons are of special interest. Each event is also ac- 
companied by at least one tau. 

We further investigated numerically single smuon pro- 
duction for X2jk 7^ within two representative fi LSP 
scenarios, i.e. for two sets of B3 mSUGRA parame- 
ters. We include the 2-body fi LSP decays via the 
generated A233 couplings in our analysis. The cross 
sections for like-sign dimuon final states are given in 
Tab. |V] and Tab. IVII and those for final states with 
three or four muons in Tab. IVII] and Tab. lVIlTl For ex- 
ample, we found resulting cross sections for exclusive 
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like-sign dimuon events of 0(100 fb) for Aj^-^^Igut = 
0.01. Additional three- and four-muon events can oc- 
cur with the same rate. This is a novel discovery mech- 
anism for the LHC and should be investigated in more 
detail, also by the LHC experimental groups. 
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APPENDIX A: CROSS SECTIONS AND 
BRANCHING RATIOS RELEVANT FOR 
SLEPTON PRODUCTION AND DECAY 

In this Appendix we give the necessary cross sections 
and branching ratios to calculate rates of all possible 
decay signatures for single slepton production at the 
LHC, within the B3 sets A and B with a fi LSP, cf. 

Eq. mm . 

In Tables IIXI and |Xl all hadronic production cross 
sections of resonant single sleptons within parameter 
Set A and Set B, respectively, are given. We consider 
here A- = 0.01 |gut, but the cross section scales with 



I A' 



ijk I 



The running of A^^^. is taken into account ac- 
cording to Eq. pi.3ip . leading to the following values 
at the SUSY scale Qsusy, cf. Eq. (lIOTll : 



Set A: A', 



2]k 



23k 



2j3 



A' 



233 



Set B: A!, 



2]k 



^2 

A' 



23k 



2j3 



233 



0.0282, 
0.0258, 
0.0281, 
0.0255, 

0.0274, 
0.0249, 
0.0269, 
0.0238, 



A' 



3jk 



A' 



33fc 



A3i3 - 
A' 



333 



A' 



3jk 



33k 



3j3 



A' 



333 



0.0282, 
0.0257, 
0.0280, 
0.0254; 

0.0271, 
: 0.0247, 
0.0266, 
0.0236, 



(A.l) 



(A.2) 



where j,k — 1,2 and Qsusy = 
Q,^,^ = 1209 GcV for Set B. 



893 GeV for Set A and 



The production cross sections include NLO SUSY- 
QCD corrections [63] . The latter depend on the 
trilincar quark-squark-slcpton coupling, hj^k , defined 
in Rcf. Numerically, it is h^k = -23.4 GeV 

(-21.2 GcV) within Set A (Set B) at the SUSY scale. 



We incorporated the running of hj^k by using the one- 
loop contributions from gauge interactions [lO|. 

Second, for the calculation of the rate for a given signa- 
ture of resonant single slepton production, the branch- 
ing ratios for the slepton decay and for the subsequent 
decay chains down to the fi LSP are needed. For 
all dominant A^^j, couplings these branching ratios are 
universal within parameter Set A and Set B. respec- 
tively, and are given in Tab. IXII 



Finally, we show in Table IXH] (Table |XlTI| all branch- 
ing ratios of fi LSP decays for different couplings Aj^j. 
at the GUT scale. Branching ratios within scenarios 
with X'ljj^ 7^ are analogous and can be obtained from 
the tables by replacing ^ by e in the final state signa- 
tures. 

In the case of a non- vanishing Ag^-^,, the fi LSP directly 
couples to the dominant L^QjDk operator and decays 
predominantly via the inverse production process, see 
also the discussion in Sect. IIII Al For the special case 
of Aggj. and mf^ < mt, however, the fi decays 
into a W boson and two jets, cf. Eq. pil.45|) . The 
corresponding matrix element and partial width are 
calculated in Appendix iBl 



APPENDIX B: THE B3 SLEPTON DECAY 

i- ^w-bdk 

A non-vanishing LiQ^Dk operator allows for slepton 
decay into a top quark t and a down-type quark dk of 
generation k, 



tdk 



(B.l) 



However, this decay mode is kinematically only al- 
lowed if > mt + rrid^. . For < mt + m^^. , the 
slepton decays via a virtual top quark. 



W-bdk. 



(B.2) 



This 3-body decay has not been considered in the lit- 
erature yet and is not implemented in the R-parity 
violating version of Herwig, either. We complete the 
picture by calculating the 3-body decay (|B.2[) in the 
following. 

The relevant parts of the supcrsymmctric Lagrangian 
are ^ 



^L^QsD, = K3kLl0i,p dkPLt + h.c. 



c 



bWt 



where Lq^ is the slepton mixing matrix, a the 
left /right eigenstate, and /3 the mass eigenstate. From 
Eq. (jB.3p . the squared matrix element (summed over 
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•^prod. 


[IDJ 








^ prod. 


[IDJ 


















^2 







^1 


~* 






— U.Ui|GUT 


z / UU 


104U 


loDU 


loDU 


zDzU 


loUU 


/I Q/1 




1 on 


1 on 


\' 


— U.Ui|GUT 


zoo 


ZDo 


/1 1 n 




ZDUU 


ZDUU 


D4.0 


D4.0 


/101 

4zi 


4zi 


\' 

^il2 


— U.Ui|GUT 


zioU 






DUz 


zuyu 


40i 


Odd 


iUo 


i4DU 


DiD 


^i21 


= O.OllcuT 


405 


1050 


602 


1430 


393 


1020 


91.9 


197 


616 


1460 


Kli 


= O.OIIgut 


1240 


220 


788 


292 


1210 


214 


216 


51.3 


806 


299 


-^'i23 


= O.OIIgut 


119 


119 


191 


191 


116 


116 


30.0 


30.0 


196 


196 




= O.OIIgut 






247 


666 










253 


681 


•^i32 


= O.OIIgut 






161 


161 










166 


166 


•^'i33 


= O.OIIgut 






69.3 


69.3 










71.1 


71.1 



TABLE IX: Complete list of hadronic cross sections for resonant single slepton/sneutrino production via A^j^j = 0.01|gut 
at the pp collider LHC (\/5 = 14 TeV) within the parameter Set A. The cross sections include QCD and SUSY-QCD 
corrections at NLO [63 . For A^3fe , sleptons cannot be produced because of the vanishing top-quark density in the proton. 





Set 




f^prod. 


[fb] 








^ prod. 


[fb] 








B 






Ki^ 






^2" 












= O.OIIgut 


885 


476 


559 


559 


949 


515 


1168 


750 


657 


657 


•^i22 


= O.OIIgut 


67.3 


67.3 


102 


102 


74.7 


74.7 


192 


192 


124 


124 


•^'il2 


= O.OIIgut 


681 


123 


414 


155 


735 


136 


976 


301 


490 


187 




= O.OIIgut 


105 


309 


155 


414 


117 


337 


269 


548 


187 


490 




= O.OIIgut 


370 


54.6 


214 


70.2 


401 


60.6 


572 


146 


255 


85.4 


^i2'i 


= O.OIIgut 


28.2 


28.2 


44.4 


44.4 


31.4 


31.4 


87.2 


87.2 


54.3 


54.3 


^VA\ 


= O.OIIgut 






60.4 


184 










73.5 


219 


-^i32 


= O.OIIgut 






38.2 


38.2 










46.7 


46.7 


-^^33 


= O.OIIgut 






14.8 


14.8 










18.2 


18.2 



TABLE X: Same as Tab.|lX]but for parameter Set B. 



final state polarizations and colors) can be derived, 
2 _ 3 A;2,L2^g2 



quark. 




Myv-hf 



mi-mfy+4:(W-b) + 



'(B.4) 

We denote the particle four-momenta by the particle 
letter, and rrif, m?,, and mw, arc the top, bottom and 
W mass, respectively. Tt is the total width of the top 



From the squared matrix element (jB.4[) we obtain 
easily the partial width for the 3-body decay (|B.2p . 
see e.g. (Qa]. We show in Fig. [12] the partial width 
T{eL — > W~bd) as a function of the left-handed selec- 
tron mass mg^ . Here we take X'l^^i = 0.01 and Ln = 1, 
in Eq. (jRl)) . 



In comparison to the 3-body decay (jB.2[) . the possible 
4-body decays via A^g^ are negligible. For example 
for the parameter Set B with non- vanishing X'^^i, the 
branching ratio of the 3-body fi LSP decay (|B.2p is 
larger by five orders of magnitude than the branching 
ratio of the 4-body fi LSP decays. 
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+ Uj dk 
dj dk 



X? 
X? 
X? 



-R ' 



' 1 



^2" 






^ f r ft" 








Uj dk 






Vt 






—> dj dk 



BRs [%] 



le-06 



X'2,k = O.OIIgut 
Set A Set B 



91.1 
8.9 
91.7 
9.3 



91.3 
8.7 
91.5 
8.4 



36.0 
7.0 
7.0 



45.7 
2.2 
2.1 



54.3 
45.7 



64.1 
35.9 



58.4 
22.5 
19.1 
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37.8 



14.7 
41.8 
43.5 

13.6 
86.4 



Set A 



O.OIIgut 
Set B 



100 
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100 



100 



36.0 
7.0 
7.0 



45.7 
2.2 
2.1 



54.3 
45.7 



64.1 
35.9 



55.5 
21.4 
18.1 
5.0 
58.8 
35.8 
5.4 



14.5 
41.2 
42.9 

1.3 
13.4 
85.2 

1.4 



TABLE XI: Table of branching ratios, BRs, that are rel- 
evant for single slepton production and decays within the 
B3 mSUGRA scenarios Set A and Set B. Two different 
non-zero B3 coupUngs are considered, X'2jk = 0.01|gut for 



columns 2 and 3 and X'^jk = 0.01|gut for columns 4 and 
5. The branching ratios for X[ji^ 7^ can be obtained from 
those for X'2jk 7^ by interchanging muon and electron fla- 
vor in the flrst four decay channels. The branching ratios 
for e_L (i^e, Br) in scenarios with X'^jk 7^ 0, i 7^ 1 are equal 
to those of (11 {i'ti, Ah) with X'^j^ ^ 0. The branching ra- 
tios for fi LSP decays are listed separately in Tabs. IXIll 
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^1 
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(2.7%) 
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TABLE XII: Branching ratios of the fi LSP for different non-zero X'^j^ couplings at the GUT scale. The branching 
ratios are calculated within the mSUGRA parameter Set A for the SUSY breaking scale Qsusy = 893 GeV. We assume 
down-type (up-type) quark mixing. Branching ratios for non-vanishing X'ljf. are analogous, with jj. replaced by e. 



Set 
B 


^1 - 


A 


-,_ A 
^1 ^ 




A' _ _ 7 

T fj, Ujdk 


A' 

^1 


T fJ,^Uj dk 


f 1 ^ r U;_idjdk 


A' 

Tl 


T v^djdk 


-^211 


49.0% 


(48.6%) 


1.7% 


(1.7%) 


(0.1%) 


0.1% 


(0.4%) 


(0.1%) 


0.1% 


(0.5%) 


-^212 


49.1% 


(-) 


1.7% 


(-) 


(5.6%) 




(41.1%) 


(6.3%) 




(46.9%) 


^213 


49.0% 


(-) 


1.7% 


(-) 


(5.7%) 


0.1% 


(41.0%) 


(6.4%) 


0.1% 


(46.9%) 


-^221 


49.1% 


(-) 


1.7% 


(-) 


(5.6%) 




(41.0%) 


(6.3%) 




(47.0%) 


^^222 


49.1% 


(49.1%) 


1.7% 


(1.7%) 


(-) 




(-) 


(-) 




(-) 


^223 


49.1% 


{') 


1.7% 


(-) 


(5.7%) 




(41.0%) 


(6.4%) 




(47.0%) 


-'^231 


49.1% 


(-) 


1.7% 


(-) 


(-) 




(-) 


(12.0%) 


0.1% 


(88.0%) 


^232 


49.1% 


(-) 


1.7% 


(-) 


(-) 




(-) 


(12.0%) 




(88.0%) 


-^233 


49.1% 


(49.1%) 


1.7% 


(1.7%) 


(-) 




(-) 


(-) 




(-) 



TABLE XIII; Branching ratios of the fi LSP for different non-zero X'2jk couplings at the GUT scale. The branching 
ratios are calculated within the mSUGRA parameter Set B for the SUSY breaking scale Qsusy = 1209 GeV. We assume 
down-type (up-type) quark mixing. Branching ratios for non-vanishing A^j, are analogous, with /i replaced by e. 
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